PANCHANGE_descriptives
Katie Thompson, Topher Huebel, Molly Davies, Kirstin Purves
17 December, 2020
This is the analysis script to prepare the data to run analyses for the project “Anxiety, depression and trauma symptom change during the COVID-19 pandemic: retrospective versus objective assessment” - Young et al (2020)
Script written by K Purves, K Thompson, C Huebel and M Davies. Email: kirstin.purves@kcl.ac.uk, katie.thompson@kcl.ac.uk, christopher.1.huebel@kcl.ac.uk, molly.davies@kcl.ac.uk
#Set up
Clear global environment
Retrieve the current date to use it for file endings to not overwrite files when one exports files
#Packages Install packages (if they are not available in your version of R)
Load packages
Registered S3 method overwritten by 'pryr':
method from
print.bytes RcppFor best results, restart R session and update pander using devtools:: or remotes::install_github('rapporter/pander')
Attaching package: 'polycor'The following object is masked from 'package:psych':
polyserialcorrplot 0.84 loaded── Attaching packages ─────────────────────────────────────── tidyverse 1.3.0 ──✓ ggplot2 3.3.2 ✓ purrr 0.3.4
✓ tibble 3.0.4 ✓ dplyr 1.0.2
✓ tidyr 1.1.2 ✓ stringr 1.4.0
✓ readr 1.3.1 ✓ forcats 0.5.0── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
x ggplot2::%+%() masks psych::%+%()
x ggplot2::alpha() masks psych::alpha()
x dplyr::filter() masks stats::filter()
x dplyr::lag() masks stats::lag()
x tibble::view() masks summarytools::view()Colour palettes
Define colours for plotting this are the standard coping colours
Choose in this chunk which palette to use 04.07.2020 - Default to 2 colour COPING palette
ggplot theme
Set up ggplot theme for the plots
Data import
Source data file paths
Read in data
[1] 65223 643Important vectors with collections of columns
Vectors for sum_score
[1] "gad.sum_score_prepan" "pcl.sum_score_prepan"
[3] "phq.sum_score_prepan" "phq.sum_score_8items_prepan"
[5] "ocir.sum_score_prepan" "gad.sum_score_base"
[7] "gad.sum_score_retro" "pcl.sum_score_base"
[9] "pcl.sum_score_retro" "phq.sum_score_base"
[11] "phq.sum_score_8items_base" "phq.sum_score_retro"
[13] "phq.sum_score_8items_retro" "ocir.sum_score_base"
[15] "ocir.sum_score_retro" [1] "phq.sum_score_prepan" "phq.sum_score_base" "phq.sum_score_retro" [1] "gad.sum_score_prepan" "gad.sum_score_base" "gad.sum_score_retro" [1] "ocir.sum_score_prepan" "ocir.sum_score_base" "ocir.sum_score_retro" [1] "pcl.sum_score_prepan" "pcl.sum_score_base" "pcl.sum_score_retro" Vectors for diff_score
[1] "phq.diff_score_retro_prepan" "phq.diff_score_base_prepan"
[3] "phq.diff_score_base_retro" "gad.diff_score_retro_prepan"
[5] "gad.diff_score_base_prepan" "gad.diff_score_base_retro"
[7] "pcl.diff_score_base_prepan" "ocir.diff_score_retro_prepan"
[9] "ocir.diff_score_base_prepan" "ocir.diff_score_base_retro" [1] "phq.diff_score_retro_prepan" "phq.diff_score_base_prepan"
[3] "phq.diff_score_base_retro" [1] "gad.diff_score_retro_prepan" "gad.diff_score_base_prepan"
[3] "gad.diff_score_base_retro" [1] "ocir.diff_score_retro_prepan" "ocir.diff_score_base_prepan"
[3] "ocir.diff_score_base_retro" Exclusions
Exclude GLAD/EDGI individuals that didn’t take part in COPING Exclude individuals with NA for the MHD categories Individuals with NAs for age/gender have already been excluded in pre-processing script, but they get excluded here to be safe and to double check previous programming.
[1] 65223 643[1] 34465 643[1] 30758frequencies and descriptives in exclusion data
pandemic change gad
Frequencies
dat$gad.pandemic_felt_feelings
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------- ------- --------- -------------- --------- --------------
Much worse 5720 16.62 16.62 16.60 16.60
A little worse 12774 37.12 53.75 37.06 53.66
No different 13640 39.64 93.39 39.58 93.24
A little better 1558 4.53 97.91 4.52 97.76
Much better 718 2.09 100.00 2.08 99.84
<NA> 55 0.16 100.00
Total 34465 100.00 100.00 100.00 100.00pandemic change phq
Frequencies
dat$phq.pandemic_felt_feelings
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------- ------- --------- -------------- --------- --------------
Much worse 5481 15.98 15.98 15.90 15.90
A little worse 13745 40.07 56.04 39.88 55.78
No different 12470 36.35 92.39 36.18 91.97
A little better 1799 5.24 97.64 5.22 97.19
Much better 811 2.36 100.00 2.35 99.54
<NA> 159 0.46 100.00
Total 34465 100.00 100.00 100.00 100.00Ethnicity
Frequencies
dat$Ethnicity_collapsed
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------- ------- --------- -------------- --------- --------------
European 30201 95.82 95.82 87.63 87.63
Mixed or multiple ethnic origins 571 1.81 97.63 1.66 89.28
Asian or Asian British 394 1.25 98.88 1.14 90.43
African or African British 159 0.50 99.39 0.46 90.89
Other 193 0.61 100.00 0.56 91.45
<NA> 2947 8.55 100.00
Total 34465 100.00 100.00 100.00 100.00Gender
Frequencies
dat$Gender_collapsed
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------------- ------- --------- -------------- --------- --------------
Male 9698 28.14 28.14 28.14 28.14
Female 24387 70.76 98.90 70.76 98.90
Non-binary/Self-defined 380 1.10 100.00 1.10 100.00
<NA> 0 0.00 100.00
Total 34465 100.00 100.00 100.00 100.00age categorical
Frequencies
dat$age_category_collapsed
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
16-18 1055 3.06 3.06 3.06 3.06
19-25 2617 7.59 10.65 7.59 10.65
26-35 5112 14.83 25.49 14.83 25.49
36-45 4781 13.87 39.36 13.87 39.36
46-55 6988 20.28 59.63 20.28 59.63
56-65 8015 23.26 82.89 23.26 82.89
66-70 3189 9.25 92.14 9.25 92.14
71-75 2010 5.83 97.97 5.83 97.97
76+ 698 2.03 100.00 2.03 100.00
<NA> 0 0.00 100.00
Total 34465 100.00 100.00 100.00 100.00Employment
Frequencies
dat$employment_prior_covid
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------- ------- --------- -------------- --------- --------------
Employed 8923 27.14 27.14 25.89 25.89
Key worker 12559 38.20 65.34 36.44 62.33
Retired 8034 24.44 89.78 23.31 85.64
Student 1507 4.58 94.36 4.37 90.01
Unemployed 1854 5.64 100.00 5.38 95.39
<NA> 1588 4.61 100.00
Total 34465 100.00 100.00 100.00 100.00anxiety disorders
Frequencies
dat$anxiety_disorders
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------- ------- --------- -------------- --------- --------------
No anxiety disorder 17177 52.97 52.97 49.84 49.84
Anxiety disorder 15251 47.03 100.00 44.25 94.09
<NA> 2037 5.91 100.00
Total 34465 100.00 100.00 100.00 100.00Depressive disorders
Frequencies
dat$depressive_disorders
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------------------- ------- --------- -------------- --------- --------------
No depressive disorder 15261 45.77 45.77 44.28 44.28
Depressive disorder 18084 54.23 100.00 52.47 96.75
<NA> 1120 3.25 100.00
Total 34465 100.00 100.00 100.00 100.00Depression and anxiety
Frequencies
dat$depression_and_anxiety
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------------------- ------- --------- -------------- --------- --------------
No depressive or anxiety disorder 12892 39.76 39.76 37.41 37.41
Depressive and anxiety disorder 12963 39.97 79.73 37.61 75.02
Only depressive disorder 4285 13.21 92.94 12.43 87.45
Only anxiety disorder 2288 7.06 100.00 6.64 94.09
<NA> 2037 5.91 100.00
Total 34465 100.00 100.00 100.00 100.00Eating disorders
Frequencies
dat$eating_disorders
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------ ------- --------- -------------- --------- --------------
No eating disorder 30748 93.54 93.54 89.22 89.22
Eating disorder 2122 6.46 100.00 6.16 95.37
<NA> 1595 4.63 100.00
Total 34465 100.00 100.00 100.00 100.00OCD
Frequencies
dat$obsessive_compulsive_disorders
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------- ------- --------- -------------- --------- --------------
No obsessive compulsive disorder 30948 92.81 92.81 89.80 89.80
Obsessive compulsive disorder 2397 7.19 100.00 6.95 96.75
<NA> 1120 3.25 100.00
Total 34465 100.00 100.00 100.00 100.00Psychosis
Frequencies
dat$psychotic_disorders
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------- ------- --------- -------------- --------- --------------
No psychotic disorder 32239 98.20 98.20 93.54 93.54
Psychotic disorder 590 1.80 100.00 1.71 95.25
<NA> 1636 4.75 100.00
Total 34465 100.00 100.00 100.00 100.00Bipolar
Frequencies
dat$mhd.mania_hypomania_bipolar_or_manicdepression
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------------------------------------- ------- --------- -------------- --------- --------------
Not Mania, hypomania, bipolar or manic-depression 32156 96.43 96.43 93.30 93.30
Mania, hypomania, bipolar or manic-depression 1189 3.57 100.00 3.45 96.75
<NA> 1120 3.25 100.00
Total 34465 100.00 100.00 100.00 100.00PTSD
Frequencies
dat$mhd.posttraumatic_stress_disorder_ptsd
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------------------------------- ------- --------- -------------- --------- --------------
Not Post-traumatic stress disorder (PTSD) 30556 91.64 91.64 88.66 88.66
Post-traumatic stress disorder (PTSD) 2789 8.36 100.00 8.09 96.75
<NA> 1120 3.25 100.00
Total 34465 100.00 100.00 100.00 100.00ASD
Frequencies
dat$autism_spectrum_disorder
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------------- ------- --------- -------------- --------- --------------
No autism spectrum disorder 32202 98.09 98.09 93.43 93.43
Autism spectrum disorder 627 1.91 100.00 1.82 95.25
<NA> 1636 4.75 100.00
Total 34465 100.00 100.00 100.00 100.00ADHD
Frequencies
dat$mhd.attention_deficit_hyperactivity_disorder
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------------------------------------------------------- ------- --------- -------------- --------- --------------
Not Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD) 32480 98.94 98.94 94.24 94.24
Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD) 349 1.06 100.00 1.01 95.25
<NA> 1636 4.75 100.00
Total 34465 100.00 100.00 100.00 100.00Personality disorders
Frequencies
dat$mhd.personality_disorder
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------------ ------- --------- -------------- --------- --------------
Not Personality disorder 31501 95.95 95.95 91.40 91.40
Personality disorder 1328 4.05 100.00 3.85 95.25
<NA> 1636 4.75 100.00
Total 34465 100.00 100.00 100.00 100.00Controls
Frequencies
dat$control
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------- ------- --------- -------------- --------- --------------
No control 21348 64.16 64.16 61.94 61.94
Control 11925 35.84 100.00 34.60 96.54
<NA> 1192 3.46 100.00
Total 34465 100.00 100.00 100.00 100.00Dont know or prefer not to say
Frequencies
dat$DidNotRespondMHD
Type: Factor
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------- ------- --------- -------------- --------- --------------
Prefer not to say 1086 3.15 3.15 3.15 3.15
Responded to MHD 33379 96.85 100.00 96.85 100.00
<NA> 0 0.00 100.00
Total 34465 100.00 100.00 100.00 100.00Data availability per questionnaire GAD, PHQ, PCL, and OCI-R
++MD: These ‘data availability’ chunks caused my R to crash, so I’m hashing them out for now. Others can test if it works for them and unhash as needed, otherwise I think they should be removed.
Data availability GAD
Data availability PHQ
Data availability PCL
Data availability OCI-R
Whole sample
Participants per Sample
Eating Disorders Genetics Initative Genetic Links to Anxiety and Depression National BioResource Rapid Assessment of MP
Frequencies
dat$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
EDGI 65 0.19 0.19 0.19 0.19
GLAD 12653 36.71 36.90 36.71 36.90
NBR 13096 38.00 74.90 38.00 74.90
RAMP 8651 25.10 100.00 25.10 100.00
<NA> 0 0.00 100.00
Total 34465 100.00 100.00 100.00 100.00Psychiatric disorders
Frequencies
dat$disorders_total_count
Type: Numeric
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
0 12386 38.85 38.85 35.94 35.94
1 5527 17.33 56.18 16.04 51.97
2 8340 26.16 82.34 24.20 76.17
3 3337 10.47 92.81 9.68 85.86
4 1393 4.37 97.17 4.04 89.90
5 569 1.78 98.96 1.65 91.55
6 239 0.75 99.71 0.69 92.24
7 64 0.20 99.91 0.19 92.43
8 23 0.07 99.98 0.07 92.49
9 6 0.02 100.00 0.02 92.51
<NA> 2581 7.49 100.00
Total 34465 100.00 100.00 100.00 100.00#HR
#SB - Impute no change prepandemic scores
[1] 21 11 16 10 1 4 13 0 14 11 0 9 16 8 7 14 14 6 4 5 9 9 21 19 4
[26] 6 7 9 7 14 [1] NA 2 9 5 NA 21 NA NA NA 3 NA 3 NA 3 3 13 NA NA 3 0 NA 5 15 17 NA
[26] NA 2 NA 4 4 [1] No different Much worse A little worse A little worse No different
[6] Much better No different No different No different A little worse
5 Levels: Much worse A little worse No different ... Much better [1] 3 1 2 2 3 5 3 3 3 2 [1] 21 11 16 10 1 4 13 0 14 11 0 9 16 8 7 14 14 6 4 5 9 9 21 19 4
[26] 6 7 9 7 14 [1] 21 2 9 5 1 21 13 0 14 3 0 3 16 3 3 13 14 6 3 0 9 5 15 17 4
[26] 6 2 9 4 4 [1] 22 9 15 11 0 9 23 5 23 11 1 20 18 15 4 11 15 19 12 10 17 11 24 27 6
[26] 9 3 11 3 8 [1] 23 1 7 7 NA 20 18 NA 23 5 NA 3 NA 7 3 11 NA NA 12 2 NA 6 19 27 NA
[26] NA 1 NA NA 3 [1] A little worse Much worse A little worse A little worse
[5] No different A little better A little worse No different
[9] A little worse A little worse
5 Levels: Much worse A little worse No different ... Much better [1] 2 1 2 2 3 4 2 3 2 2 [1] 22 9 15 11 0 9 23 5 23 11 1 20 18 15 4 11 15 19 12 10 17 11 24 27 6
[26] 9 3 11 3 8 [1] 23 1 7 7 0 20 18 5 23 5 1 3 18 7 3 11 15 19 12 2 17 6 19 27 6
[26] 9 1 11 3 3#Baseline survey completion (N per month)
[1] "2020-04-30 20:51:40 UTC" "2020-04-30 18:55:12 UTC"
[3] "2020-04-30 19:02:04 UTC" "2020-04-30 18:59:37 UTC"
[5] "2020-05-01 10:02:34 UTC" "2020-04-30 22:16:01 UTC"
[7] "2020-04-30 20:36:37 UTC" "2020-04-30 19:02:03 UTC"
[9] "2020-04-30 18:57:33 UTC" "2020-05-01 08:54:01 UTC"
[11] "2020-04-30 21:47:03 UTC" "2020-04-30 19:10:42 UTC"
[13] "2020-04-30 19:29:36 UTC" "2020-04-30 19:20:42 UTC"
[15] "2020-05-01 15:09:23 UTC" "2020-04-30 18:59:55 UTC"
[17] "2020-05-04 16:40:21 UTC" "2020-05-01 16:06:45 UTC"
[19] "2020-05-19 10:34:16 UTC" "2020-04-30 19:11:26 UTC"Vector
Correlation matrix
Warning in polychor(x, y, ML = ML, std.err = std.err): inadmissible correlation
set to -0.9999
Warning in polychor(x, y, ML = ML, std.err = std.err): inadmissible correlation
set to -0.9999Warning in hetcor.data.frame(as.data.frame(dat[, mhd.items]), use =
"pairwise.complete.obs"): the correlation matrix has been adjusted to make it
positive-definite
Two-Step Estimates
Correlations/Type of Correlation:
depressive_disorders
depressive_disorders 1
anxiety_disorders 0.8377
eating_disorders 0.5343
obsessive_compulsive_disorders 0.5269
psychotic_disorders 0.5033
mhd.mania_hypomania_bipolar_or_manicdepression 0.5367
mhd.posttraumatic_stress_disorder_ptsd 0.5807
autism_spectrum_disorder 0.3882
mhd.attention_deficit_hyperactivity_disorder 0.3696
mhd.personality_disorder 0.6383
control -0.8096
anxiety_disorders
depressive_disorders Polychoric
anxiety_disorders 1
eating_disorders 0.4969
obsessive_compulsive_disorders 0.6313
psychotic_disorders 0.4568
mhd.mania_hypomania_bipolar_or_manicdepression 0.431
mhd.posttraumatic_stress_disorder_ptsd 0.5833
autism_spectrum_disorder 0.4385
mhd.attention_deficit_hyperactivity_disorder 0.4147
mhd.personality_disorder 0.5838
control -0.8075
eating_disorders
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders 1
obsessive_compulsive_disorders 0.5707
psychotic_disorders 0.3641
mhd.mania_hypomania_bipolar_or_manicdepression 0.3653
mhd.posttraumatic_stress_disorder_ptsd 0.3867
autism_spectrum_disorder 0.2962
mhd.attention_deficit_hyperactivity_disorder 0.3309
mhd.personality_disorder 0.5191
control -0.6858
obsessive_compulsive_disorders
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders 1
psychotic_disorders 0.3917
mhd.mania_hypomania_bipolar_or_manicdepression 0.3917
mhd.posttraumatic_stress_disorder_ptsd 0.4481
autism_spectrum_disorder 0.3976
mhd.attention_deficit_hyperactivity_disorder 0.3908
mhd.personality_disorder 0.502
control -0.7346
psychotic_disorders
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders 1
mhd.mania_hypomania_bipolar_or_manicdepression 0.7086
mhd.posttraumatic_stress_disorder_ptsd 0.4476
autism_spectrum_disorder 0.2901
mhd.attention_deficit_hyperactivity_disorder 0.3375
mhd.personality_disorder 0.6044
control -0.701
mhd.mania_hypomania_bipolar_or_manicdepression
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression 1
mhd.posttraumatic_stress_disorder_ptsd 0.3492
autism_spectrum_disorder 0.2813
mhd.attention_deficit_hyperactivity_disorder 0.342
mhd.personality_disorder 0.5855
control -0.6882
mhd.posttraumatic_stress_disorder_ptsd
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd 1
autism_spectrum_disorder 0.2916
mhd.attention_deficit_hyperactivity_disorder 0.2851
mhd.personality_disorder 0.5077
control -0.6867
autism_spectrum_disorder
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd Polychoric
autism_spectrum_disorder 1
mhd.attention_deficit_hyperactivity_disorder 0.5726
mhd.personality_disorder 0.3777
control -0.6313
mhd.attention_deficit_hyperactivity_disorder
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd Polychoric
autism_spectrum_disorder Polychoric
mhd.attention_deficit_hyperactivity_disorder 1
mhd.personality_disorder 0.3698
control -0.6354
mhd.personality_disorder
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd Polychoric
autism_spectrum_disorder Polychoric
mhd.attention_deficit_hyperactivity_disorder Polychoric
mhd.personality_disorder 1
control -0.7725
control
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd Polychoric
autism_spectrum_disorder Polychoric
mhd.attention_deficit_hyperactivity_disorder Polychoric
mhd.personality_disorder Polychoric
control 1
Standard Errors/Numbers of Observations:
depressive_disorders
depressive_disorders 33345
anxiety_disorders 0.004037
eating_disorders 0.01226
obsessive_compulsive_disorders 0.0117
psychotic_disorders 0.02142
mhd.mania_hypomania_bipolar_or_manicdepression 0.01569
mhd.posttraumatic_stress_disorder_ptsd 0.01078
autism_spectrum_disorder 0.02049
mhd.attention_deficit_hyperactivity_disorder 0.02606
mhd.personality_disorder 0.01498
control Inf
anxiety_disorders
depressive_disorders 32428
anxiety_disorders 32428
eating_disorders 0.01209
obsessive_compulsive_disorders 0.01054
psychotic_disorders 0.0202
mhd.mania_hypomania_bipolar_or_manicdepression 0.01546
mhd.posttraumatic_stress_disorder_ptsd 0.01032
autism_spectrum_disorder 0.01985
mhd.attention_deficit_hyperactivity_disorder 0.02525
mhd.personality_disorder 0.01397
control Inf
eating_disorders
depressive_disorders 32836
anxiety_disorders 31919
eating_disorders 32870
obsessive_compulsive_disorders 0.01299
psychotic_disorders 0.02372
mhd.mania_hypomania_bipolar_or_manicdepression 0.01942
mhd.posttraumatic_stress_disorder_ptsd 0.01538
autism_spectrum_disorder 0.02513
mhd.attention_deficit_hyperactivity_disorder 0.02924
mhd.personality_disorder 0.0157
control 1.851
obsessive_compulsive_disorders
depressive_disorders 33345
anxiety_disorders 32428
eating_disorders 32836
obsessive_compulsive_disorders 33345
psychotic_disorders 0.02248
mhd.mania_hypomania_bipolar_or_manicdepression 0.01825
mhd.posttraumatic_stress_disorder_ptsd 0.01397
autism_spectrum_disorder 0.02216
mhd.attention_deficit_hyperactivity_disorder 0.0268
mhd.personality_disorder 0.01562
control 11.73
psychotic_disorders
depressive_disorders 32804
anxiety_disorders 31887
eating_disorders 32826
obsessive_compulsive_disorders 32804
psychotic_disorders 32829
mhd.mania_hypomania_bipolar_or_manicdepression 0.01621
mhd.posttraumatic_stress_disorder_ptsd 0.02076
autism_spectrum_disorder 0.03517
mhd.attention_deficit_hyperactivity_disorder 0.03917
mhd.personality_disorder 0.01887
control 3.285
mhd.mania_hypomania_bipolar_or_manicdepression
depressive_disorders 33345
anxiety_disorders 32428
eating_disorders 32836
obsessive_compulsive_disorders 33345
psychotic_disorders 32804
mhd.mania_hypomania_bipolar_or_manicdepression 33345
mhd.posttraumatic_stress_disorder_ptsd 0.01847
autism_spectrum_disorder 0.02949
mhd.attention_deficit_hyperactivity_disorder 0.03287
mhd.personality_disorder 0.01632
control 1.964
mhd.posttraumatic_stress_disorder_ptsd
depressive_disorders 33345
anxiety_disorders 32428
eating_disorders 32836
obsessive_compulsive_disorders 33345
psychotic_disorders 32804
mhd.mania_hypomania_bipolar_or_manicdepression 33345
mhd.posttraumatic_stress_disorder_ptsd 33345
autism_spectrum_disorder 0.02388
mhd.attention_deficit_hyperactivity_disorder 0.02903
mhd.personality_disorder 0.01514
control 1.983
autism_spectrum_disorder
depressive_disorders 32804
anxiety_disorders 31887
eating_disorders 32826
obsessive_compulsive_disorders 32804
psychotic_disorders 32829
mhd.mania_hypomania_bipolar_or_manicdepression 32804
mhd.posttraumatic_stress_disorder_ptsd 32804
autism_spectrum_disorder 32829
mhd.attention_deficit_hyperactivity_disorder 0.02798
mhd.personality_disorder 0.02533
control 3.353
mhd.attention_deficit_hyperactivity_disorder
depressive_disorders 32804
anxiety_disorders 31887
eating_disorders 32826
obsessive_compulsive_disorders 32804
psychotic_disorders 32829
mhd.mania_hypomania_bipolar_or_manicdepression 32804
mhd.posttraumatic_stress_disorder_ptsd 32804
autism_spectrum_disorder 32829
mhd.attention_deficit_hyperactivity_disorder 32829
mhd.personality_disorder 0.03046
control 2.038
mhd.personality_disorder control
depressive_disorders 32804 33271
anxiety_disorders 31887 32382
eating_disorders 32826 32800
obsessive_compulsive_disorders 32804 33271
psychotic_disorders 32829 32777
mhd.mania_hypomania_bipolar_or_manicdepression 32804 33271
mhd.posttraumatic_stress_disorder_ptsd 32804 33271
autism_spectrum_disorder 32829 32777
mhd.attention_deficit_hyperactivity_disorder 32829 32777
mhd.personality_disorder 32829 32777
control 1.749 33273Create the correlation matrix plot using the corrplot package 
#Correlations of retro impute and prepan ##Create function to use for all of the corplots
###Whole sample prepan retro correlation
Frequencies
prepan_retro_dat$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
EDGI 65 0.51 0.51 0.51 0.51
GLAD 12653 99.49 100.00 99.49 100.00
<NA> 0 0.00 100.00
Total 12718 100.00 100.00 100.00 100.00
Two-Step Estimates
Correlations/Type of Correlation:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 1 Pearson
gad.sum_score_retro 0.5143 1
phq.sum_score_prepan 0.6904 0.4485
phq.sum_score_retro 0.4443 0.7445
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan Pearson Pearson
gad.sum_score_retro Pearson Pearson
phq.sum_score_prepan 1 Pearson
phq.sum_score_retro 0.5886 1
Standard Errors/Numbers of Observations:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 12692 9063
gad.sum_score_retro 0.007727 9082
phq.sum_score_prepan 0.004652 0.008394
phq.sum_score_retro 0.008274 0.004961
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan 12662 9411
gad.sum_score_retro 9058 8074
phq.sum_score_prepan 12681 9405
phq.sum_score_retro 0.006739 9431
P-values for Tests of Bivariate Normality:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan
gad.sum_score_retro 6.754e-309
phq.sum_score_prepan 2.342e-223 2.296e-275
phq.sum_score_retro 4.022e-131 7.718e-221
phq.sum_score_prepan
gad.sum_score_prepan
gad.sum_score_retro
phq.sum_score_prepan
phq.sum_score_retro 1.101e-125Warning in corCi(x = x, keys = keys, n.iter = n.iter, p = p, overlap = overlap, :
Finding confidence intervals from a correlation matrix, but n is not specified, arbitrarily set to 1000
Frequencies
whole_retro_prepan_scatter$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
EDGI 65 0.51 0.51 0.51 0.51
GLAD 12653 99.49 100.00 99.49 100.00
<NA> 0 0.00 100.00
Total 12718 100.00 100.00 100.00 100.00Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 33 rows containing non-finite values (stat_smooth).Warning: Removed 33 rows containing missing values (geom_point).
Warning: Removed 33 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 42 rows containing non-finite values (stat_smooth).Warning: Removed 42 rows containing missing values (geom_point).
Warning: Removed 42 rows containing missing values (geom_point).
###Little/much worse prepan retro correlation
[1] "gad.sum_score_prepan" "gad.sum_score_retro" "phq.sum_score_prepan"
[4] "phq.sum_score_retro" Frequencies
prepan_retro_dat_worse$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------ --------- -------------- --------- --------------
EDGI 42 0.54 0.54 0.54 0.54
GLAD 7795 99.46 100.00 99.46 100.00
<NA> 0 0.00 100.00
Total 7837 100.00 100.00 100.00 100.00
Two-Step Estimates
Correlations/Type of Correlation:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 1 Pearson
gad.sum_score_retro 0.5476 1
phq.sum_score_prepan 0.69 0.4785
phq.sum_score_retro 0.4592 0.7366
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan Pearson Pearson
gad.sum_score_retro Pearson Pearson
phq.sum_score_prepan 1 Pearson
phq.sum_score_retro 0.5976 1
Standard Errors/Numbers of Observations:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 7820 7808
gad.sum_score_retro 0.007924 7825
phq.sum_score_prepan 0.005931 0.008729
phq.sum_score_retro 0.009446 0.005473
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan 7805 6982
gad.sum_score_retro 7805 6988
phq.sum_score_prepan 7817 6979
phq.sum_score_retro 0.007697 6997
P-values for Tests of Bivariate Normality:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan
gad.sum_score_retro 8.019e-142
phq.sum_score_prepan 1.022e-144 2.245e-111
phq.sum_score_retro 1.966e-123 1.629e-103
phq.sum_score_prepan
gad.sum_score_prepan
gad.sum_score_retro
phq.sum_score_prepan
phq.sum_score_retro 2.61e-134Warning in corCi(x = x, keys = keys, n.iter = n.iter, p = p, overlap = overlap, :
Finding confidence intervals from a correlation matrix, but n is not specified, arbitrarily set to 1000
Frequencies
little_much_worse_retro_prepan$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
EDGI 65 0.51 0.51 0.51 0.51
GLAD 12653 99.49 100.00 99.49 100.00
<NA> 0 0.00 100.00
Total 12718 100.00 100.00 100.00 100.00Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 14 rows containing non-finite values (stat_smooth).Warning: Removed 14 rows containing missing values (geom_point).
Warning: Removed 14 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 15 rows containing non-finite values (stat_smooth).Warning: Removed 15 rows containing missing values (geom_point).
Warning: Removed 15 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 29 rows containing non-finite values (stat_smooth).Warning: Removed 29 rows containing missing values (geom_point).
Warning: Removed 29 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 11 rows containing non-finite values (stat_smooth).Warning: Removed 11 rows containing missing values (geom_point).
Warning: Removed 11 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 21 rows containing non-finite values (stat_smooth).Warning: Removed 21 rows containing missing values (geom_point).
Warning: Removed 21 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 32 rows containing non-finite values (stat_smooth).Warning: Removed 32 rows containing missing values (geom_point).
Warning: Removed 32 rows containing missing values (geom_point).
###Little/much better prepan retro correlation
[1] "gad.sum_score_prepan" "gad.sum_score_retro" "phq.sum_score_prepan"
[4] "phq.sum_score_retro" Frequencies
prepan_retro_dat_better$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------ --------- -------------- --------- --------------
EDGI 7 0.56 0.56 0.56 0.56
GLAD 1252 99.44 100.00 99.44 100.00
<NA> 0 0.00 100.00
Total 1259 100.00 100.00 100.00 100.00
Two-Step Estimates
Correlations/Type of Correlation:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 1 Pearson
gad.sum_score_retro 0.4858 1
phq.sum_score_prepan 0.6512 0.4084
phq.sum_score_retro 0.3567 0.6977
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan Pearson Pearson
gad.sum_score_retro Pearson Pearson
phq.sum_score_prepan 1 Pearson
phq.sum_score_retro 0.5346 1
Standard Errors/Numbers of Observations:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 1257 1255
gad.sum_score_retro 0.02158 1257
phq.sum_score_prepan 0.01628 0.02355
phq.sum_score_retro 0.02651 0.01559
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan 1253 1085
gad.sum_score_retro 1253 1086
phq.sum_score_prepan 1255 1084
phq.sum_score_retro 0.02171 1087
P-values for Tests of Bivariate Normality:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan
gad.sum_score_retro 2.716e-29
phq.sum_score_prepan 7.975e-34 3.25e-08
phq.sum_score_retro 1.024e-27 1.832e-07
phq.sum_score_prepan
gad.sum_score_prepan
gad.sum_score_retro
phq.sum_score_prepan
phq.sum_score_retro 1.275e-12Warning in corCi(x = x, keys = keys, n.iter = n.iter, p = p, overlap = overlap, :
Finding confidence intervals from a correlation matrix, but n is not specified, arbitrarily set to 1000
Frequencies
little_much_better_retro_prepan$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
EDGI 65 0.51 0.51 0.51 0.51
GLAD 12653 99.49 100.00 99.49 100.00
<NA> 0 0.00 100.00
Total 12718 100.00 100.00 100.00 100.00Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 3 rows containing non-finite values (stat_smooth).Warning: Removed 3 rows containing missing values (geom_point).
Warning: Removed 3 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 1 rows containing non-finite values (stat_smooth).Warning: Removed 1 rows containing missing values (geom_point).
Warning: Removed 1 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 4 rows containing non-finite values (stat_smooth).Warning: Removed 4 rows containing missing values (geom_point).
Warning: Removed 4 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 4 rows containing non-finite values (stat_smooth).
Warning: Removed 4 rows containing missing values (geom_point).
Warning: Removed 4 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 6 rows containing non-finite values (stat_smooth).Warning: Removed 6 rows containing missing values (geom_point).
Warning: Removed 6 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h`geom_smooth()` using formula 'y ~ x'Warning: Removed 10 rows containing non-finite values (stat_smooth).Warning: Removed 10 rows containing missing values (geom_point).
Warning: Removed 10 rows containing missing values (geom_point).
Create long data frame
Disorders per Gender
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
100% bar per diagnosis with percentages of contribution per sample
Disorders per sample
| Disorder_inc | Sample | n | Prop |
|---|---|---|---|
| Anxiety disorders | EDGI | 38 | 0.00 |
| Anxiety disorders | GLAD | 10098 | 0.66 |
| Anxiety disorders | NBR | 1765 | 0.12 |
| Anxiety disorders | RAMP | 3350 | 0.22 |
| Depressive disorders | EDGI | 47 | 0.00 |
| Depressive disorders | GLAD | 11455 | 0.63 |
| Depressive disorders | NBR | 2756 | 0.15 |
| Depressive disorders | RAMP | 3826 | 0.21 |
| Eating disorders | EDGI | 61 | 0.03 |
| Eating disorders | GLAD | 1384 | 0.65 |
| Eating disorders | NBR | 195 | 0.09 |
| Eating disorders | RAMP | 482 | 0.23 |
| OCDs | EDGI | 18 | 0.01 |
| OCDs | GLAD | 1627 | 0.68 |
| OCDs | NBR | 145 | 0.06 |
| OCDs | RAMP | 607 | 0.25 |
| Psychotic disorders | GLAD | 401 | 0.68 |
| Psychotic disorders | NBR | 37 | 0.06 |
| Psychotic disorders | RAMP | 152 | 0.26 |
| Bipolar disorders | EDGI | 3 | 0.00 |
| Bipolar disorders | GLAD | 874 | 0.74 |
| Bipolar disorders | NBR | 72 | 0.06 |
| Bipolar disorders | RAMP | 240 | 0.20 |
| PTSD | EDGI | 6 | 0.00 |
| PTSD | GLAD | 1851 | 0.66 |
| PTSD | NBR | 231 | 0.08 |
| PTSD | RAMP | 701 | 0.25 |
| Autism | EDGI | 6 | 0.01 |
| Autism | GLAD | 384 | 0.61 |
| Autism | NBR | 46 | 0.07 |
| Autism | RAMP | 191 | 0.30 |
| ADHD | EDGI | 2 | 0.01 |
| ADHD | GLAD | 201 | 0.58 |
| ADHD | NBR | 32 | 0.09 |
| ADHD | RAMP | 114 | 0.33 |
| Personality disorders | EDGI | 6 | 0.00 |
| Personality disorders | GLAD | 951 | 0.72 |
| Personality disorders | NBR | 76 | 0.06 |
| Personality disorders | RAMP | 295 | 0.22 |
| Controls | EDGI | 1 | 0.00 |
| Controls | GLAD | 139 | 0.01 |
| Controls | NBR | 9247 | 0.78 |
| Controls | RAMP | 2538 | 0.21 |
Warning: Ignoring unknown parameters: binwidth
Disorder per ethnicity
Warning: Ignoring unknown parameters: binwidth
Age categories
Age bar plot for whole sample (not grouped) Summary dataset for age and display it as a bar plot
| age_category | n | Prop | order |
|---|---|---|---|
| 16-18 | 1055 | 0.03 | 1 |
| 19-25 | 2617 | 0.08 | 2 |
| 26-35 | 5112 | 0.15 | 3 |
| 36-45 | 4781 | 0.14 | 4 |
| 46-55 | 6988 | 0.20 | 5 |
| 56-65 | 8015 | 0.23 | 6 |
| 66-70 | 3189 | 0.09 | 7 |
| 71-75 | 2010 | 0.06 | 8 |
| 76-80 | 544 | 0.02 | 9 |
| 81-85 | 116 | 0.00 | 10 |
| 86-90 | 27 | 0.00 | 11 |
| 91-100 | 6 | 0.00 | 12 |
| 100+ | 5 | 0.00 | 13 |
Education
Pie charts
summary dataset with frequencies
| highest_education | n | freq |
|---|---|---|
| GCSE/CSE | 4224 | 0.12 |
| NVQ | 2148 | 0.06 |
| A-levels | 6742 | 0.20 |
| University | 19948 | 0.58 |
pie chart labels
Warning: Ignoring unknown parameters: face
Warning: Ignoring unknown parameters: face
Warning: Ignoring unknown aesthetics: face
pie chart - no labels
Ethnicity_collapsed
Pie charts
| Ethnicity_collapsed | n | freq |
|---|---|---|
| European | 30201 | 0.96 |
| Mixed or multiple ethnic origins | 571 | 0.02 |
| Asian or Asian British | 394 | 0.01 |
| African or African British | 159 | 0.01 |
| Other | 193 | 0.01 |
pie chart with only European label
Warning: Ignoring unknown parameters: face
Pie chart - no labels
Gender collapsed
Pie charts
| Gender_collapsed | n | freq |
|---|---|---|
| Male | 9698 | 0.28 |
| Female | 24387 | 0.71 |
| Non-binary/Self-defined | 380 | 0.01 |
Warning: Ignoring unknown parameters: face
Pie chart - only labelling two largest groups (male and female)
Warning: Ignoring unknown parameters: face
Pie chart - no labels 
Compare the GLAD/EDGI and COPING subcohorts: EDGI, GLAD, NBR, RAMP
+++CH: Diagnosis, ethnicity, gender, age
age_category
Proportion tables
| age_category | Sample | n | Prop |
|---|---|---|---|
| 16-18 | EDGI | 7 | 0.00 |
| 16-18 | GLAD | 341 | 0.01 |
| 16-18 | NBR | 10 | 0.00 |
| 16-18 | RAMP | 697 | 0.02 |
| 19-25 | EDGI | 17 | 0.00 |
| 19-25 | GLAD | 1537 | 0.04 |
| 19-25 | NBR | 291 | 0.01 |
| 19-25 | RAMP | 772 | 0.02 |
| 26-35 | EDGI | 21 | 0.00 |
| 26-35 | GLAD | 2850 | 0.08 |
| 26-35 | NBR | 1165 | 0.03 |
| 26-35 | RAMP | 1076 | 0.03 |
| 36-45 | EDGI | 8 | 0.00 |
| 36-45 | GLAD | 2459 | 0.07 |
| 36-45 | NBR | 1440 | 0.04 |
| 36-45 | RAMP | 874 | 0.03 |
| 46-55 | EDGI | 4 | 0.00 |
| 46-55 | GLAD | 2821 | 0.08 |
| 46-55 | NBR | 2644 | 0.08 |
| 46-55 | RAMP | 1519 | 0.04 |
| 56-65 | EDGI | 5 | 0.00 |
| 56-65 | GLAD | 1908 | 0.06 |
| 56-65 | NBR | 3901 | 0.11 |
| 56-65 | RAMP | 2201 | 0.06 |
| 66-70 | EDGI | 3 | 0.00 |
| 66-70 | GLAD | 433 | 0.01 |
| 66-70 | NBR | 1967 | 0.06 |
| 66-70 | RAMP | 786 | 0.02 |
| 71-75 | GLAD | 231 | 0.01 |
| 71-75 | NBR | 1265 | 0.04 |
| 71-75 | RAMP | 514 | 0.01 |
| 76-80 | GLAD | 59 | 0.00 |
| 76-80 | NBR | 334 | 0.01 |
| 76-80 | RAMP | 151 | 0.00 |
| 81-85 | GLAD | 13 | 0.00 |
| 81-85 | NBR | 61 | 0.00 |
| 81-85 | RAMP | 42 | 0.00 |
| 86-90 | GLAD | 1 | 0.00 |
| 86-90 | NBR | 13 | 0.00 |
| 86-90 | RAMP | 13 | 0.00 |
| 91-100 | NBR | 2 | 0.00 |
| 91-100 | RAMP | 4 | 0.00 |
| 100+ | NBR | 3 | 0.00 |
| 100+ | RAMP | 2 | 0.00 |
| Sample | age_category | n | Prop |
|---|---|---|---|
| EDGI | 16-18 | 7 | 0.11 |
| EDGI | 19-25 | 17 | 0.26 |
| EDGI | 26-35 | 21 | 0.32 |
| EDGI | 36-45 | 8 | 0.12 |
| EDGI | 46-55 | 4 | 0.06 |
| EDGI | 56-65 | 5 | 0.08 |
| EDGI | 66-70 | 3 | 0.05 |
| GLAD | 16-18 | 341 | 0.03 |
| GLAD | 19-25 | 1537 | 0.12 |
| GLAD | 26-35 | 2850 | 0.23 |
| GLAD | 36-45 | 2459 | 0.19 |
| GLAD | 46-55 | 2821 | 0.22 |
| GLAD | 56-65 | 1908 | 0.15 |
| GLAD | 66-70 | 433 | 0.03 |
| GLAD | 71-75 | 231 | 0.02 |
| GLAD | 76-80 | 59 | 0.00 |
| GLAD | 81-85 | 13 | 0.00 |
| GLAD | 86-90 | 1 | 0.00 |
| NBR | 16-18 | 10 | 0.00 |
| NBR | 19-25 | 291 | 0.02 |
| NBR | 26-35 | 1165 | 0.09 |
| NBR | 36-45 | 1440 | 0.11 |
| NBR | 46-55 | 2644 | 0.20 |
| NBR | 56-65 | 3901 | 0.30 |
| NBR | 66-70 | 1967 | 0.15 |
| NBR | 71-75 | 1265 | 0.10 |
| NBR | 76-80 | 334 | 0.03 |
| NBR | 81-85 | 61 | 0.00 |
| NBR | 86-90 | 13 | 0.00 |
| NBR | 91-100 | 2 | 0.00 |
| NBR | 100+ | 3 | 0.00 |
| RAMP | 16-18 | 697 | 0.08 |
| RAMP | 19-25 | 772 | 0.09 |
| RAMP | 26-35 | 1076 | 0.12 |
| RAMP | 36-45 | 874 | 0.10 |
| RAMP | 46-55 | 1519 | 0.18 |
| RAMP | 56-65 | 2201 | 0.25 |
| RAMP | 66-70 | 786 | 0.09 |
| RAMP | 71-75 | 514 | 0.06 |
| RAMP | 76-80 | 151 | 0.02 |
| RAMP | 81-85 | 42 | 0.00 |
| RAMP | 86-90 | 13 | 0.00 |
| RAMP | 91-100 | 4 | 0.00 |
| RAMP | 100+ | 2 | 0.00 |
Cross tables
| age_category | n_EDGI | n_GLAD | n_NBR | n_RAMP | Prop_EDGI | Prop_GLAD | Prop_NBR | Prop_RAMP |
|---|---|---|---|---|---|---|---|---|
| 16-18 | 7 | 341 | 10 | 697 | 0.11 | 0.03 | 0.00 | 0.08 |
| 19-25 | 17 | 1537 | 291 | 772 | 0.26 | 0.12 | 0.02 | 0.09 |
| 26-35 | 21 | 2850 | 1165 | 1076 | 0.32 | 0.23 | 0.09 | 0.12 |
| 36-45 | 8 | 2459 | 1440 | 874 | 0.12 | 0.19 | 0.11 | 0.10 |
| 46-55 | 4 | 2821 | 2644 | 1519 | 0.06 | 0.22 | 0.20 | 0.18 |
| 56-65 | 5 | 1908 | 3901 | 2201 | 0.08 | 0.15 | 0.30 | 0.25 |
| 66-70 | 3 | 433 | 1967 | 786 | 0.05 | 0.03 | 0.15 | 0.09 |
| 71-75 | NA | 231 | 1265 | 514 | NA | 0.02 | 0.10 | 0.06 |
| 76-80 | NA | 59 | 334 | 151 | NA | 0.00 | 0.03 | 0.02 |
| 81-85 | NA | 13 | 61 | 42 | NA | 0.00 | 0.00 | 0.00 |
| 86-90 | NA | 1 | 13 | 13 | NA | 0.00 | 0.00 | 0.00 |
| 91-100 | NA | NA | 2 | 4 | NA | NA | 0.00 | 0.00 |
| 100+ | NA | NA | 3 | 2 | NA | NA | 0.00 | 0.00 |
| age_category | EDGI | GLAD | NBR | RAMP |
|---|---|---|---|---|
| 16-18 | 7 | 341 | 10 | 697 |
| 19-25 | 17 | 1537 | 291 | 772 |
| 26-35 | 21 | 2850 | 1165 | 1076 |
| 36-45 | 8 | 2459 | 1440 | 874 |
| 46-55 | 4 | 2821 | 2644 | 1519 |
| 56-65 | 5 | 1908 | 3901 | 2201 |
| 66-70 | 3 | 433 | 1967 | 786 |
| 71-75 | 0 | 231 | 1265 | 514 |
| 76-80 | 0 | 59 | 334 | 151 |
| 81-85 | 0 | 13 | 61 | 42 |
| 86-90 | 0 | 1 | 13 | 13 |
| 91-100 | 0 | 0 | 2 | 4 |
| 100+ | 0 | 0 | 3 | 2 |
Chi square test
Warning in chisq.test(.): Chi-squared approximation may be incorrect| statistic | p.value | parameter | method |
|---|---|---|---|
| 5856.536 | 0 | 36 | Pearson’s Chi-squared test |
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
Ethnicity_collapsed
Proportion tables
| Ethnicity_collapsed | Sample | n | Prop |
|---|---|---|---|
| European | EDGI | 60 | 0.00 |
| European | GLAD | 12090 | 0.35 |
| European | NBR | 9927 | 0.29 |
| European | RAMP | 8124 | 0.24 |
| Mixed or multiple ethnic origins | EDGI | 3 | 0.00 |
| Mixed or multiple ethnic origins | GLAD | 276 | 0.01 |
| Mixed or multiple ethnic origins | NBR | 118 | 0.00 |
| Mixed or multiple ethnic origins | RAMP | 174 | 0.01 |
| Asian or Asian British | EDGI | 1 | 0.00 |
| Asian or Asian British | GLAD | 98 | 0.00 |
| Asian or Asian British | NBR | 126 | 0.00 |
| Asian or Asian British | RAMP | 169 | 0.00 |
| African or African British | GLAD | 37 | 0.00 |
| African or African British | NBR | 46 | 0.00 |
| African or African British | RAMP | 76 | 0.00 |
| Other | EDGI | 1 | 0.00 |
| Other | GLAD | 122 | 0.00 |
| Other | RAMP | 70 | 0.00 |
| NA | GLAD | 30 | 0.00 |
| NA | NBR | 2879 | 0.08 |
| NA | RAMP | 38 | 0.00 |
| Sample | Ethnicity_collapsed | n | Prop |
|---|---|---|---|
| EDGI | European | 60 | 0.92 |
| EDGI | Mixed or multiple ethnic origins | 3 | 0.05 |
| EDGI | Asian or Asian British | 1 | 0.02 |
| EDGI | Other | 1 | 0.02 |
| GLAD | European | 12090 | 0.96 |
| GLAD | Mixed or multiple ethnic origins | 276 | 0.02 |
| GLAD | Asian or Asian British | 98 | 0.01 |
| GLAD | African or African British | 37 | 0.00 |
| GLAD | Other | 122 | 0.01 |
| GLAD | NA | 30 | 0.00 |
| NBR | European | 9927 | 0.76 |
| NBR | Mixed or multiple ethnic origins | 118 | 0.01 |
| NBR | Asian or Asian British | 126 | 0.01 |
| NBR | African or African British | 46 | 0.00 |
| NBR | NA | 2879 | 0.22 |
| RAMP | European | 8124 | 0.94 |
| RAMP | Mixed or multiple ethnic origins | 174 | 0.02 |
| RAMP | Asian or Asian British | 169 | 0.02 |
| RAMP | African or African British | 76 | 0.01 |
| RAMP | Other | 70 | 0.01 |
| RAMP | NA | 38 | 0.00 |
Cross tables
| Ethnicity_collapsed | n_EDGI | n_GLAD | n_NBR | n_RAMP | Prop_EDGI | Prop_GLAD | Prop_NBR | Prop_RAMP |
|---|---|---|---|---|---|---|---|---|
| European | 60 | 12090 | 9927 | 8124 | 0.92 | 0.96 | 0.76 | 0.94 |
| Mixed or multiple ethnic origins | 3 | 276 | 118 | 174 | 0.05 | 0.02 | 0.01 | 0.02 |
| Asian or Asian British | 1 | 98 | 126 | 169 | 0.02 | 0.01 | 0.01 | 0.02 |
| Other | 1 | 122 | NA | 70 | 0.02 | 0.01 | NA | 0.01 |
| African or African British | NA | 37 | 46 | 76 | NA | 0.00 | 0.00 | 0.01 |
| NA | NA | 30 | 2879 | 38 | NA | 0.00 | 0.22 | 0.00 |
| Ethnicity_collapsed | EDGI | GLAD | NBR | RAMP |
|---|---|---|---|---|
| European | 60 | 12090 | 9927 | 8124 |
| Mixed or multiple ethnic origins | 3 | 276 | 118 | 174 |
| Asian or Asian British | 1 | 98 | 126 | 169 |
| African or African British | 0 | 37 | 46 | 76 |
| Other | 1 | 122 | 0 | 70 |
| NA | 0 | 30 | 2879 | 38 |
Chi square test
Warning in chisq.test(.): Chi-squared approximation may be incorrect| statistic | p.value | parameter | method |
|---|---|---|---|
| 5106.32 | 0 | 15 | Pearson’s Chi-squared test |
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
Gender_collapsed
Proportion tables
| Gender_collapsed | Sample | n | Prop |
|---|---|---|---|
| Male | EDGI | 1 | 0.00 |
| Male | GLAD | 2279 | 0.07 |
| Male | NBR | 5641 | 0.16 |
| Male | RAMP | 1777 | 0.05 |
| Female | EDGI | 62 | 0.00 |
| Female | GLAD | 10127 | 0.29 |
| Female | NBR | 7416 | 0.22 |
| Female | RAMP | 6782 | 0.20 |
| Non-binary/Self-defined | EDGI | 2 | 0.00 |
| Non-binary/Self-defined | GLAD | 247 | 0.01 |
| Non-binary/Self-defined | NBR | 39 | 0.00 |
| Non-binary/Self-defined | RAMP | 92 | 0.00 |
| Sample | Gender_collapsed | n | Prop |
|---|---|---|---|
| EDGI | Male | 1 | 0.02 |
| EDGI | Female | 62 | 0.95 |
| EDGI | Non-binary/Self-defined | 2 | 0.03 |
| GLAD | Male | 2279 | 0.18 |
| GLAD | Female | 10127 | 0.80 |
| GLAD | Non-binary/Self-defined | 247 | 0.02 |
| NBR | Male | 5641 | 0.43 |
| NBR | Female | 7416 | 0.57 |
| NBR | Non-binary/Self-defined | 39 | 0.00 |
| RAMP | Male | 1777 | 0.21 |
| RAMP | Female | 6782 | 0.78 |
| RAMP | Non-binary/Self-defined | 92 | 0.01 |
Cross tables
| Gender_collapsed | n_EDGI | n_GLAD | n_NBR | n_RAMP | Prop_EDGI | Prop_GLAD | Prop_NBR | Prop_RAMP |
|---|---|---|---|---|---|---|---|---|
| Male | 1 | 2279 | 5641 | 1777 | 0.02 | 0.18 | 0.43 | 0.21 |
| Female | 62 | 10127 | 7416 | 6782 | 0.95 | 0.80 | 0.57 | 0.78 |
| Non-binary/Self-defined | 2 | 247 | 39 | 92 | 0.03 | 0.02 | 0.00 | 0.01 |
| Gender_collapsed | EDGI | GLAD | NBR | RAMP |
|---|---|---|---|---|
| Male | 1 | 2279 | 5641 | 1777 |
| Female | 62 | 10127 | 7416 | 6782 |
| Non-binary/Self-defined | 2 | 247 | 39 | 92 |
Chi square test
Warning in chisq.test(.): Chi-squared approximation may be incorrect| statistic | p.value | parameter | method |
|---|---|---|---|
| 2455.701 | 0 | 6 | Pearson’s Chi-squared test |
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
highest_education
Proportion tables
| highest_education | Sample | n | Prop |
|---|---|---|---|
| GCSE/CSE | EDGI | 7 | 0.00 |
| GCSE/CSE | GLAD | 1078 | 0.03 |
| GCSE/CSE | NBR | 1876 | 0.06 |
| GCSE/CSE | RAMP | 1263 | 0.04 |
| NVQ | EDGI | 3 | 0.00 |
| NVQ | GLAD | 918 | 0.03 |
| NVQ | NBR | 1227 | 0.04 |
| A-levels | EDGI | 20 | 0.00 |
| A-levels | GLAD | 2713 | 0.08 |
| A-levels | NBR | 2320 | 0.07 |
| A-levels | RAMP | 1689 | 0.05 |
| University | EDGI | 35 | 0.00 |
| University | GLAD | 7515 | 0.23 |
| University | NBR | 6983 | 0.21 |
| University | RAMP | 5415 | 0.16 |
| Sample | highest_education | n | Prop |
|---|---|---|---|
| EDGI | GCSE/CSE | 7 | 0.11 |
| EDGI | NVQ | 3 | 0.05 |
| EDGI | A-levels | 20 | 0.31 |
| EDGI | University | 35 | 0.54 |
| GLAD | GCSE/CSE | 1078 | 0.09 |
| GLAD | NVQ | 918 | 0.08 |
| GLAD | A-levels | 2713 | 0.22 |
| GLAD | University | 7515 | 0.61 |
| NBR | GCSE/CSE | 1876 | 0.15 |
| NBR | NVQ | 1227 | 0.10 |
| NBR | A-levels | 2320 | 0.19 |
| NBR | University | 6983 | 0.56 |
| RAMP | GCSE/CSE | 1263 | 0.15 |
| RAMP | A-levels | 1689 | 0.20 |
| RAMP | University | 5415 | 0.65 |
Cross tables
| highest_education | n_EDGI | n_GLAD | n_NBR | n_RAMP | Prop_EDGI | Prop_GLAD | Prop_NBR | Prop_RAMP |
|---|---|---|---|---|---|---|---|---|
| GCSE/CSE | 7 | 1078 | 1876 | 1263 | 0.11 | 0.09 | 0.15 | 0.15 |
| NVQ | 3 | 918 | 1227 | NA | 0.05 | 0.08 | 0.10 | NA |
| A-levels | 20 | 2713 | 2320 | 1689 | 0.31 | 0.22 | 0.19 | 0.20 |
| University | 35 | 7515 | 6983 | 5415 | 0.54 | 0.61 | 0.56 | 0.65 |
| highest_education | EDGI | GLAD | NBR | RAMP |
|---|---|---|---|---|
| GCSE/CSE | 7 | 1078 | 1876 | 1263 |
| NVQ | 3 | 918 | 1227 | 0 |
| A-levels | 20 | 2713 | 2320 | 1689 |
| University | 35 | 7515 | 6983 | 5415 |
Chi square test
Warning in chisq.test(.): Chi-squared approximation may be incorrect| statistic | p.value | parameter | method |
|---|---|---|---|
| 1207.855 | 0 | 12 | Pearson’s Chi-squared test |
Bar plots
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
key_worker
Proportion tables
| key_worker | Sample | n | Prop |
|---|---|---|---|
| No key worker | EDGI | 42 | 0.00 |
| No key worker | GLAD | 7249 | 0.21 |
| No key worker | NBR | 8660 | 0.25 |
| No key worker | RAMP | 5841 | 0.17 |
| Key worker | EDGI | 23 | 0.00 |
| Key worker | GLAD | 5366 | 0.16 |
| Key worker | NBR | 4414 | 0.13 |
| Key worker | RAMP | 2765 | 0.08 |
| NA | GLAD | 38 | 0.00 |
| NA | NBR | 22 | 0.00 |
| NA | RAMP | 45 | 0.00 |
| Sample | key_worker | n | Prop |
|---|---|---|---|
| EDGI | No key worker | 42 | 0.65 |
| EDGI | Key worker | 23 | 0.35 |
| GLAD | No key worker | 7249 | 0.57 |
| GLAD | Key worker | 5366 | 0.42 |
| GLAD | NA | 38 | 0.00 |
| NBR | No key worker | 8660 | 0.66 |
| NBR | Key worker | 4414 | 0.34 |
| NBR | NA | 22 | 0.00 |
| RAMP | No key worker | 5841 | 0.68 |
| RAMP | Key worker | 2765 | 0.32 |
| RAMP | NA | 45 | 0.01 |
Cross tables
| key_worker | n_EDGI | n_GLAD | n_NBR | n_RAMP | Prop_EDGI | Prop_GLAD | Prop_NBR | Prop_RAMP |
|---|---|---|---|---|---|---|---|---|
| No key worker | 42 | 7249 | 8660 | 5841 | 0.65 | 0.57 | 0.66 | 0.68 |
| Key worker | 23 | 5366 | 4414 | 2765 | 0.35 | 0.42 | 0.34 | 0.32 |
| NA | NA | 38 | 22 | 45 | NA | 0.00 | 0.00 | 0.01 |
| key_worker | EDGI | GLAD | NBR | RAMP |
|---|---|---|---|---|
| No key worker | 42 | 7249 | 8660 | 5841 |
| Key worker | 23 | 5366 | 4414 | 2765 |
| NA | 0 | 38 | 22 | 45 |
Chi square test
Warning in chisq.test(.): Chi-squared approximation may be incorrect| statistic | p.value | parameter | method |
|---|---|---|---|
| 332.6547 | 0 | 6 | Pearson’s Chi-squared test |
Bar plots
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth